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Electromagnetic interaction of a magnetized rotating star 
with a conducting disk 
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Abstract A conducting disk significantly changes the 
generation of the electromagnetic radiation excited by 
the rotation of the magnetic field frozen to a star. Due 
to the reflection of waves from a disk there appear waves 
propagating toward a star, not only outward a star as 
it takes place for the magneto-dipole radiation. Be- 
cause that the angular momentum can be transformed 
from a disk to a star when the inner edge of a disk 
approaches the light surface of a rotating star. This 
is purely electromagnetic effect. At some distance of a 
disk from a star, = r* ~ c/u) 8 , the stellar angular 
momentum losses due to the electromagnetic radiation 
become zero. It results the stable stellar rotation. 
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1 Introduction 

The interaction of an accretion disk with a magnetized 
star is very complicated and is not understandable yet. 
The observations of X-ray binaries demonstrate many 
physical phenomenon [1]. It seems useful to solve a set 
of simple physical problems throwing light on the na- 
ture of interaction of an accretion flow with a star. One 
of these problems is the interaction of a rotating mag- 
netic field with a conducting matter of a disk. The phe- 
nomena depends on the conductivity of the disk matter 
a. The physics is different for the high and low conduc- 
tivities. Aly [2] consider the disk conductivity is high 
and the magnetic field does not penetrate to a disk. 
Opposite, Bardou and Heyvaerts [3] consider the con- 
ductivity is low due to abnormal resistivity of a turbu- 
lent disk plasma, and the field penetrate to a disk freely. 
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The interesting combination of penetration and repul- 
sion of magnetic field lines was suggested by Lovelace 
at. al [4]. We consider here a disk to be stationary, 
without accretion and with the classical properties. 

The ionized plasma of an accretion disk has the con- 
ductivity a = 10 13 (T e /lel/) 3 / 2 (A/I0)- 1 s- 1 , which is 
high enough to consider a disk as an ideal conduc- 
tor. Here T e is the temperature of electrons in a disk 
which is larger than lOeF, A is the Coulomb loga- 
rithm. At such conductivity a the width of the skin 
layer \ s k = (rc 2 /^) 1 / 2 is less than the disk width H. 
The value of r is the characteristic time of turbulent 
motion in the a- disk, r = H/ Vk , Vk is the Kcplerian ve- 
locity of the disk rotation. The condition a ^> c 2 /Hvk 
is well fulfilled in the inner parts of a disk. 

Here we consider the oblique magnetic field rotating 
with the frequency of the star rotation u) s , and a disk 
having the infinite conductivity and finite dimensions 
. The environment of this system is the vacuum. Wc 
will see that these conditions though the boundary con- 
dition on a disk surface results to the repulsion of the 
magnetic field from a disk. 



2 Statement of problem 

We are solving the following problem: a rotating neu- 
tron star with the angular velocity w s and the magnetic 
momentum fi is surrounded by a disk (sec figure 1). The 
angle between lo s and \i is x- The matter of a disk is 
ideal, i.e. the conductivity a = oo. It means E' = 0, 
where E' is the electric field in the frame of a disk. The 
field in the laboratory frame is 



E=--[v d B]. 

c 



(1) 



Here is the rotational velocity of the disk matter 
around the a star v<j = u^e,*. We consider a disk is 
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rotating with the Keplerian velocity oc r *' 2 in the 
6 direction 



/Gil/ S \ 1/2 . 



(2) 



Here M s is the stellar mass, G is the gravitational con- 
stant and r, <j> are the spherical coordinates. To create 
the electric field (1) the matter of a disk must polarize 
and produces the charge density 
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The rotating charge density produces the toroidal elec- 
tric current = v^gsiui), 
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We consider that a disk occupies the region r c i < r < 
oo, tt/2—5 < ■& < ir/2+5, where is inner edge of a disk 
and 25 is its width in the i? direction. Integrating the 
Maxwell equation (cwZB)^ = (47r/c)(j0 i^<t> = 

inside a disk) over the angle i? from ir/2 — S to tt/2 + S 
we obtain the condition connecting two values of B r on 
both sides of a disk 



Br{\-5)-B r {^ + 5)=\ 



vl cos 2 5 
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Here the current f ^ is the toroidal electric current ex- 
citing in a disk by the infinitely small electric field E'^. 
The boundary condition (4) is the condition for the cur- 
rent j'j, arising inside a disk. Another two boundary 
conditions follow from (1) 



E^-S) = E 4 ,(^ + S) = 0,r>r d , 



(5) 




Fig. 1 A rotating neutron star with a dipole magnetic field 
is surrounded by a disk. 



(E r - — sin^)|^ =i _5 = 0, r > r d - (6) 
c 2 

(-Er " sint9B^)|^ = - +5 = 0, r > r d . 

c 2 



3 Solution 

The general solution of the wave equation in vacuum in 
the spherical coordinates consists of the sum of the mul- 
tipolc components, magnetic and electric. The mag- 
netic multipolc of the n-th order is 



B, 



Pl{cosV)z-iKJ nH +b%J_ n _i 
1 8P -i d zi\a*J n+ , +bf,J_ 



n(n+l) dti z 

Pi 1 d 



E, 



n (n+l)i&5l& Z2 \ a n J n+± + bf l J_ n _i] 
0, 

^ z" 3 [a,f J„ + 1 + 6^ J_„_ i ] , 



n(n+l) sin 
n(n+l) d-d 
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(7) 
erf,. 



Here and below P^(cos$) are the associated Legen- 
drc functions, J n +i/2( z ) is the Besscl functions of the 
( n+1/2) order, z = ru) s /c. The electric multipole of the 
n-th order is 



E„ 



&zt[a*J n+i +b%J_ n _ 



i(n+l) 9i9 z dz 
i Pi 1 d 



n(n+l) ^hilJl9l Z2 J n+i + "n ^-n- 3 J 
0, 



B, 



n(n+l) sin i9 ' 



erf,, 
(8) 

e,A. 



n(n+l) 3# 

All quantities B„ and E n are oscillating functions with 
the frequency of the stellar rotation lu s , (B,E)„ cx 
e Mt>-iu,t _ Multipoles (7,8) contain waves travelling in 
the positive direction from a neutron star to the infin- 
ity, in this case b n = i(—l) n+1 a n and the combination 

a n Jn+1/2 + b n J- n - 1/2 = a n-ff„+i/ 2 is tnc Hankel func- 
tion of the first kind 



(i) 

n+l/2 



(—) 

\2ttz J 



1/2 



Here Z n (l/z) is the polinomial of the n-th order of 
the argument Also for the wave propagating in 

the negative direction from the infinity b n = i(—l) n a n , 
and we have the Hankel function of the second kind, 



H 



(2) 



- e' lz {2/Tiz) 1 / 2 Z n {l/z). When the disk is ab- 
sent the condition at the infinity demands the relation 



n+l, 



In the presence of a disk the part 



of the wave energy can be reflected from a disk and we 
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have the mixture of both types of waves. In general 
case the coefficients a B ' E ,b B ' E are the arbitrary com- 
plex numbers and can be determined from the bound- 
ary conditions (5,6). The boundary condition on the 
stellar surface r — R s (z = z — R s uj s /c) are 

B r (z = zq) = B r (inside); E§{z = zq) = E^(inside); 

E^(z = zq) = E '^(inside). 
Fields inside are 



B = ^sm X e i 



E= ^sin X e^~ lWs * 



2P 1 1 (cosi?), 

_i_5L(l_ l z 2) 
sin 17 V 2 ' 



(9) 



-1^(1 -K), 

sintf ' 3 dti ' 
' l 3i) ' l 3 sini? 



0, e<h. 



Here is the magnetic moment of a star. We con- 
sider here only the varying component of the magnetic 
dipole, proportional to sin^, which produces an elec- 
tromagnetic radiation. From (5,6) and (9) it follows 



(of J 3/2 (z) + b? J_ 3 /a(*)) 
(a B J n+1/2 (z) + b B J- n -i/2{z)) 



O t 1 3/2 . 



d .1 



z 1,2 {a E J b/2 + b E J^ /2 ) 



(j- z Z 1,2 {a E J n+1/2 + b E J_ n _ 1/2 )j 



0,n^l; (10) 



= 2 ^3 z o sin x; 



0, n + 2. 



Introducing the notations 



^n+l/2( z ) 



g-(Z Jn+1/2) 



f(^/2j_ n _ 1/2 ) 



z=z 
(11) 

and the dimensionless quantities a n ,b n , dividing them 
on the value 2fism\/ R^, we rewrite (10) in the form 



bf = -afgf + Sl , b% = -egg, n^l, 



B R 



(12) 



Let us note that because zq <C 1 (R s -C c/w s ) the co- 
efficients gf , g E , s±, s 2 are small values. In the lowest 
order of zq they are 



T(l/2-n) 
r(3/2 + n) 



i)r(i/2 



(I) 



2)1+1 



nr(3/2- 
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Sl 



2 -3/2 r( _ 1/2)zo 3 = _2- 1 / 2 7T 1 / 2 Z 3 ; 



S2 « -2- 7 / 2 r(-3/2)z 5 = -^^V- 

Here r(x) is the Gamma function. In the case of a disk 
Eq. (12) are the connections between the coefficients 
b n and a n . When a disk is absent connections are con- 
ditions b n = i{— l)'" +1 a rl , that is the condition for the 
radiation toward the infinity when there are no reflected 
waves. Then Eq. (12) shows that nonzero coefficients 
arc af and af only 



H + gf) Sl 
(i + ef 2 ) 



(* + Q§)S2 



-ISi 



132 = ~3 [ ~ 



1/2 



1/2 



2 (1 + 9?) 
It means that a rotating magnetized star without a disk 
radiates the magnetodipole radiation with the ampli- 
tude af ~ Zq and also the electroquadrupole radia- 
tion with the amplitude a E << af . The disk changes 
this picture, and a star radiates also many multipolcs 
1,2. 

To find the coefficients af' E we need use the bound- 
ary conditions on a disk (5,6). The boundary condition 
(5) means 



Eoo 1 
n=l n(n+l) 

dP^M ( a Bj 
dtf \ a n J n+± 



Ph$o) -i d - 
. „ ' z 2 -z-z< 

sin oz 



i a n<Jn+ 



b E J 

0; z > Zd, #o = f ± 
(14) 
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We use the recurrent relations for the Bcssel functions 
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Here the values s\ and s 2 are 

z 3/2 
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Then multiplying Eq. (14) on the Gcgcnbauer polino- 
(13) mial A k ij 2 (z) and integrating it over the contour pass- 
ing through the infinity as shown on the figure 2a, we 
obtain the matrix equation connected the coefficients 
a B and a E . The Gegenbauer polinomial A k x / 2 is the 



■1 



function conjugated to the Bessel functions J n +i/2{z), figure 2b. 



<k/2 



A kti (z) = 2-Hk + -) ]T 

m=0 

and possesses the property [5] 



r(i +k-m) /z\2"»-n 



f 

2iir 



The result is 



(n+1) dtfo \Un^k,n "k,n)^ n T ( 2 „+l) sini? 

I &k + l,n &n,k-l \ I El Cfc,,i-1 Cfc,n+1 \ 

^ n n+1 )T«nl n n+1 ' 

2^siCfe,i + Hi§^S2(C fc ,i - §C fc , 3 ). 



of = 



(15) 



Here coefficients C k n are 



Cfe.n = 2~ 7 / 4,1/2 ^-n- 1/2 ^ 1/2 ^ = 
(-1)"T(± + + m) 



<k/2 



(-i) i+1 (fc + ^)E 



' m!(Z + 1 - m)!r(| + fc - m - /) 

The coefficients Cfc i7l are not equal to zero if n + k = 
21 + 1, / is the integer number. 

If we introduce two matrices K\ B n , K\ E n and the vec- 
tor a\, 



A is 

IV k,n 



a IE 
iV fc.n 



1 



■{QnCk.n — <5fc,n); 



sini? (2n+l 
l n+1 



S k -i, n 
n+1 



(16) 



then we obtain 



1 P 1 

2 fn 

10 sin 



fc,i 



rCfe,3)i 



a IB B | a IE E 1 



(17) 



The second matrix equation connecting a E and a B is 
followed from the condition (6). Doing the same proce- 
dure ( multiplying on the Gegenbauer polinomial and 
integrating in the complex plane z ), but taking into 
account that the Keplerian velocity, Vk oc z -1 / 2 , is the 
two fold function, we have to introduce the cut in the 
complex plane from the point z = Zd to the infinity. 
The resulting contour C has the shape as shown on the 



in 9P^ sini? [/ ■+ n B\l_(j+ , 

7T d$ 2n+l I \ u k,n-l ' u k,n-lUn ) n \ u k,n+l^ 

d k,n+l@n )^+T a n + P n [ 5 k,n ~ C ki nQ% + 
7T n(n+l) y u k,n a k,n^n > 

-^s 2 {Pi + fg-sj sin^o) - P%s 2 C k , 2 



E in 



dt dP} ■ a 



(18) 



Here 



4,1/2(2)4+1/2^) dz, z = z + iO; 



d fe „ = / A M/2 (z)J_„_ 1/2 (z)dz, z = z + iO, 

ft = (r g uj s /2c) 1 ' 2 << 1 is the Keplerian parameter, r s 
is the gravitational radius of a star, r g = 2GM s /c 2 , 



k = 5.5 ■ 10" 



M,\ 1/2 /P^ 1/2 



Mr 



Is 



P s is the period of a star rotation, P s = 2tt/uj s . Again, 
introducing the matrices A 2B n , A 2E n and the vector er 2 , 
we rewrite Eq. (18) in the form 



\2B B 1 \2E E „ 



(19) 



A2B i k dP,\ sin i? Q 
yl fc,n 7T dtfo 2n+l 

(^fe,n+l + d k n+1 Qn )^j" 



«n-l+ rf fe ,n-l^)^ 



A2£ = pi 



4,n ~ Cfc,„£?f + ~ n(n+1) (<„ " d fc ,nPn)J 

(20) 



or = — 



% Sl !g sin tf _& (S2 pi. 



SlTlAT sini^o 



52^^,2. 



Combining Eq. (17) and Eq. (19) we find that the 
coefficients a B and are the solution of the general 
matrix equation 

Ak, n a = 
where 7^ andl^ are the vectors 



(21) 



It 



1? 



The matrix A k>n is 



KlB klE 

A _ I Il k,n 

Ak.n I a 2B \2E 

Y1 fc,n 1Y k,n 
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4 Torque acting on a star 

Our purpose is to determine the change of the torque 
acting on a star under the influence of a rotating disk, 



K = - / r x J, x B ds, 



(22) 



where 3 S is the surface current on the star surface 



-^-x{B}, {B} 

An r 



B 



r=R s +0 



B 



r=R 3 -0 



(23) 



{B} is the discontinuous of the tangential magnetic 
field. The torque is 



K 



Rs 

c 



(24) 



where B sr is the radial magnetic field on the stellar 
surface 

Integrating in Eq. (24) over the surface we find the 
component of the torque K z along the axis of the rota- 
tion. It turns out to be proportional to two coefficients 
of the expansion of the vacuum field over multipolcs 
(7,8) 



K z = -I^Im(af)[^(z^J 3/2 ) 

i^ 1/2 J„ 3 / 2 )^ 



-V2J. 



5/2 



-3/2 

z-V*J_ 



1 s m 

Z = Z Q 

5/2£ {zl/2j ^ /2) 



(25) 

That means that only two types of radiation, the mag- 
nctodipolc and the electroquadrupole, retard the star 
rotation. The contribution of the magnetodipole com- 
ponent is much larger than the electroquadrupole one 
because Zq « 1. In the absence of a disk the torque is 



Kz = K 



AID 



2 /rsin x 3 

3 m z ° 



2 [i 2 ui 3 sin 2 x 



, (26) 



\K* ID \ = 1.5-10 30 B 2 r2 



R 



lOfcr 



sm x er 9- 



An conducting disk distorts the radiation. There ap- 
pears another multipole components of the radiation. 
The amplitudes of them are connected by the relation 
(21). The coefficients Ak, n contain also the imaginary 
parts which are proportional to the Keplerian param- 
eter k. The disk also changes the conditions at the 
infinity - a disk reflects the wave energy. As a result 
a star can as lost the angular momentum as obtain it 
from a disk. The effect depends on the disk parame- 
ter Zd = rdOJs/c. The qualitative estimation demands 



the solution of the matrix equation (21) numerically. 
The result of calculations for the usual neutron star 
parameters R = 10km, M s = 1M , ui s = Is -1 and 
Zd = 1, S = 10° is presented on the table 1. 



n 


1 


2 


3 


4 


5 


6 


Ima^ 


0.25 


0.15 


0.05 


0.013 


0.004 


0.0027 


Rea% 


0.21 


0.18 


0.05 


0.013 


0.008 


0.003 


Ima% 


0.25 


0.20 


0.075 


0.027 


0.007 


0.0027 


Rea% 


0.25 


0.28 


0.065 


0.027 


0.008 


0.0027 



Table 1. Amplitudes of n-th magnetic multipoles 
and n-th electric multipoles for the first sixth 
waves (7,8). They are normolized by their values in 
the vacuum: by Imaf = (tt /2) 1 / 2 Zq and a„ by 
/maf = -(^/72) 1 /2 Z 5. 

We see that though the main component is af , there 
appear compatible amplitudes af ,af ,af . Among elec- 
tric multipolcs the electrodipole amplitude af becomes 
general. For the larger n > 4 the amplitudes of mul- 
tipoles fall down. It permits us to restrict ourself by 
the finite dimension of vectors ~ci , 77" in equation (21), 
k = 12. The calculation of the torque K z , acting on a 
star, is presented on the figure 3. For comparison we 
draw also the result of calculation for k = 8. It is seen 
that the difference between k = 8 and k = 12 is not sig- 
nificant and our approximation to use finite numbers of 
equations in infinite matrix equation (21) is valid. The 
main result is the change of sign of K z at Zd = z* 1. 
For the chosen parameters z* = 0.8. The behavior of 
the torque K z over Zd can be well approximated the by 
the simple formula 



K 7 = K 



AID 



(27) 



This expression does not take into account the small 
maximum of K z at Zd ~ 2 — 3, when a star transmits 
its angular momentum not only to the electromagnetic 
radiation, but also to a disk. 

At Zd < z* a disk transmits the angular momen- 
tum to a star, it spins up. Let us note that z* m 1 
is just the region where the electromagnetic radiation 
forms. A disk changes radiation conditions in the zone 
where it originates. Such electromagnetic stellar spin 
up exists independently of the mechanical angular mo- 
mentum transition, which is proportional to the value 
of the accretion matter rate M. 

The structure of the magnetic field is shown on the 
figure 4. Here we present only the magnetic field pro- 
duced by the varying component of the magnetic dipolc, 
proportional to \x sin x-, which generates the electromag- 
netic radiation. 

We see that the magnetic field lines repeal from an 
ideal disk. It is the result of calculations under the 
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Fig. 2 The contour of integration in the complex plane z. 
The contour a) corresponds to the boundary condition (5), 
and the contour b) corresponds to the boundary condition 
(6) on the disk surface. 





Fig. 4 Projection of magnetic field lines of the oblique 
dipole onto the plane orthogonal to the direction of the mag- 
netic moment \x. Owing to the prjective view and the large 
scale of the order of light cylinder radius c/u a , some details 
on this figure look like peculiarities (cusps, touchs), but they 
are not in reality. 



Fig. 3 The torque K z , acting on a star, in the units of 
j^md versus mner e dge of a disk z&. 5 = 10°, a) k=12, 
b) k=8. 
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boundary conditions (5,6) on the disk surface. It is not 
postulated ad hoc as did in [2] . 

5 Discussion 

We showed that presence of a conducting disk in the 
magnctosphcrc of a rotating magnetized star changes 
significantly the electromagnetic radiation of a star if 
the inner edge of a disk n is closer than the radius of 
the stellar light cylinder c/uj s . The torque acting on a 
star K z changes its sign when < r* ~ c/u s . Instead 
of the spinning down of the star rotation, a star begins 
to spin up, getting the angular momentum from a disk. 
At the position of the inner disk radius near the light 
cylinder radius, r^ = r* , the torque K z becomes zero. 
It means that a star will not change its angular momen- 
tum and the period of its rotation will be constant. It 
can happen after the period of the gas accretion onto 
a star when it gains the angular momentum, value of 
which is proportional to the accretion rate M. As was 
shown [6], under the mass accretion onto a magnetized 
star, the torque K z has always definite sign, and there 
is no situation of the stable star rotation. When the 
accretion stops, M = 0, a star pushes out a disk from 
its vicinity, that is, so called, the propeller regime. The 
inner edge of a disk moves toward the light cylinder dis- 
tance where a disk prevents the stellar spin down due 
to the radiation of electromagnetic waves. The period 
of the stellar rotation will not change when a disk is 
presented at the distance from a star ~ c/to s . Such 
period of the star life is observed among X-rays binaries 
[!]■ 
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